Gap opening in graphene by simple periodic inhomogeneous strain 
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Using ab-initio methods, we show first that the uniform deformation either leaves graphene 
(semi) metallic or opens up a small gap yet only beyond its mechanical breaking point, contrary 
to claims in the literature based on tight-binding (TB) calculations. It is possible, however, to 
open up a global gap by a sine-like inhomogeneous deformation applied along any direction but the 
armchair one, with the largest gap along the zigzag direction (~ 1.0 eV) and without any electro- 
static gating. The gap opening has a threshold character with very sharp rise when the ratio of the 
amplitude A and the period of the sine wave deformation A exceeds (A/X) c ~ 0.1 and the inversion 
symmetry is preserved, while it is threshold-less when the symmetry is broken. The gap opening 
occurs in graphene mesh on boron-nitride substrate. 

PACS numbers: 73.22.Pr, 81.05.ue, 62.25.-g 



I. INTRODUCTION 

Graphene, a hexagonally packed single layer of carbon 
atomfl, is considered a strong candidate for post-silicon 
electronic devices. To meet this expectation, however, 
it is crucial to be able to open up a band gap at the 
Fermi level, which exists in silicon and is necessary for a 
switchable device. Since graphene has outstanding me- 
chanical properties and capable of sustaining huge atomic 
distortions 2 3 up to about 27% 4 it becomes very interest- 
ing to use strain as a tool to introduce a band gap in 
graphene. 

The idea of strain-induced band gap has been theo- 
retically explored by Pereira et a/P within the standard 
tight-binding (TB) approximation. The gap opening pro- 
cess involves merging of two inequivalent Dirac points, 
D and —Dy^* Such a process requires some ~ 23% of 
stretching along the zigzag chains 5 . The subsequent to 
Ref|5| and more realistic ab-initio calculation^"^ con- 
firmed that the strain along a zigzag direction is indeed 
the most effective in annihilating the Dirac points. How- 
ever, the ab-initio methods predicted even higher criti- 
cal tensile strain 27%) than that obtained within the 
tight-binding model. Moreover, it was discovered that 
annihilation of the Dirac points does not automatically 
lead to a gap opening, because a <r* conduction band 
quickly moves down towards the Fermi level with the 
strain. As a result, when the Dirac points merge, the en- 
ergy gap either does not ope n at a lP, or opens by a tiny 
45 meW Thus, the result ^SEHHI suggest that a tangi- 
ble gap in graphene is unlikely under uniaxial tensions 
up to the graphene failure strain of 25-27% rl 4 l n [ The 
papers^MJ considered a tension accompanied by a Pois- 
son's contraction. Using a TB approach, Cocco et alf^l 
have tried a uniform pure shear strain e xy = e yx ^ 
associated with the elastic constant C44. They claimed a 
moderate critical value of 16% and expected it to lower 
down to 12% if used along with the uniaxial strain in the 
armchair direction. 

The idea of using inhomogeneous strain fields u with 
ViUk 7^ has also been studied extensively bearing 



on TB-derived notion of a pseudo-magnetic field B cx 
y^EMU] j£ was claimed that two-dimensional corruga- 
tions with triangular symmetry open up a gap without a 
threshold^. The local gaps in graphene 'bubbles' have 
been reported in Ref. [14[ while it was speculated that the 
one-dimensional corrugations ('wrinkles') that are fre- 
quently observed in a suspended graphene could open up 
the gap but only in combination with alternating electro- 
static gating correlated with the 'wrinkles^. 

Here, we pursue two goals with the use of ab-initio cal- 
culations applied similarly to our earlier study of a flexo- 
electric [polarization induced by d x u z d y u z , (d x u z d y u z ) 2 
and (d x d y u z ) 2 ] 'sister' system B-N monolayei^. First, 
we show that contrary to Ref. 12 graphene 's band struc- 
ture remains gapless under any uniform strains not ex- 
ceeding its breaking point. Secondly, we demonstrate 
that the energy gap can be nevertheless induced mechan- 
ically by applying realistic nonuniform sine-wave defor- 
mations (similar to the 'wrinkles') in all but the armchair 
directions, and without any periodic gatinj^. The qual- 
itative difference with TB-derived continuous pseudo- 
magnetic models (see Ref. 15 and references therein) ap- 
pears to be an account for strong ppa overlap between 
neighboring p z orbitals on Carbon atoms appearing at 
any flexing with \7u ^ 0, in addition to the standard ppn 
overlap, the only one retained in TB models of graphene. 
Since the rehybridization responsible for the gap opening 
is local, the results do not depend on standard procedure 
of employing supercell for the band structure calcula- 
tions and atomic relaxation. Another feature missing in 
the continuous flexo-' magnetic' models is that they do 
not distinguish between flexing that is preserving an in- 
version symmetry versus a symmetry breaking one. 



II. HOMOGENEOUS AND ID PERIODIC 
DEFORMATIONS OF GRAPHENE 

We begin with defining the homogeneous and inho- 
mogeneous (sinusoidal) deformations of graphene. Se- 
lect the unit vectors for real lattices of an undistorted 
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2D graphene as ai = § (a/3, — 1,), CL2 — § (a/3, l) 
with the reciprocal lattice vectors G\ = (l, — a/3) , 
G2 = ^7=^ (l, a/3). They are chosen in such a way that 
the x— axis is along the armchair direction (ai + 02), 
while the y axis is along the zigzag direction (— a\ + CI2). 
Let us consider a general homogeneous deformation of a 
graphene sheet. Excluding overall translation and rota- 
tions, such a change can be described by the displace- 
ments: 



eo + V 7 
7 e - 77 



(1) 



where eo corresponds to an isotropic distortion (dilata- 
tion), while the parameters n and 7 describe the pure 
independent shear strains. By rotating the frame around 
the z axis through an angle <p given by cos (20) = 
v/ V 1 ! 2 + 7 2 7 one can diagonalize the stress tensor in |lj): 



package^ A 16x16x1 Monkhorst-Pack fe-point gricPH 
has been used in the case of flat graphene. Approxi- 
mately the same k-point density was kept in going from 
the planar to a corrugated graphene. The sheets have 
been simulated by a slab-supercell approach with the 
inter-planar distances of 30as to ensure negligible wave 
function overlap between the replica sheets. For the 
plane- wave expansion of the valence and conduction band 
wave-functions, a cutoff energy was chosen to be 80 Ry. 
The plane waves Troullier-Martins^ and Fritz-Haber- 
Institute^ pseudopotentials have been used for the cal- 
culations in the local density (LDA) and the generalized 
gradient (GGA) approximations, respectively. In both 
cases, Carbon 2s and 2p electrons have been considered 
as valence states. 



III. NO GAP IN GRAPHENE SUBJECT TO A 
HOMOGENOUS DEFORMATION 



(5x>\ (e + T \ /V\ = f eil 0\ (x'\ 

WJ V ^-r] \y>) ~ \ 622 J [J/J* 

(2) 

where r \[rf + 7 2 , en = e + r and e 2 2 = e - t. 

For inhomogeneous deformations, it is sufficient to 
consider a periodic out-of-plane atomic displacements 
u z (r) = Asin(k • r + (/?), where k is the undulation 
wave vector, r is the in-plane vector with the origin at 
the lattice inversion center, and tp some phase. To set 
up the supercell calculations, we use a periodically and 
commensurately distorted graphene sheets, the corruga- 
tion wave vector is k = 27re/A, where the unit vector e 
and the wavelength A are given by e = A/A, A(n,m) = 
na\ + ma2, A = a\fn 2 + nm + m 2 , where n and rn are 
integers, |ai| = |ot2 1 = a -> cf! 17 1 18 1 The corrugation 
with the period A(n, m) leads to a rectangular supercell 
with translational vectors A(n, m) and T = Na\ +Ma2, 
T _L A, N = (2m+n)/d R , M -(2n + m)/d R , where d R 
is the greatest common divisor of 2m + n and 2n + mP^ 
The number of graphene unit cells per supercell A x T is 
N g = 2X 2 /(a 2 d R ). 

For the undeformed graphene, the two inequivalent 
Dirac points are at (2G 1 + G 2 )/3 and (G 1 + 2G 2 )/3, 
respectively, at the two corners K and of the first 
Brillouin zone (BZ}^. It is easy to find the positions 
of these points inside the rectangular BZ of the super- 
cell (when A — >• 0) using a so-called zone-folding tech- 
nique. Namely, k = 2i\ej\ is of one of the reciprocal 
vectors of the supercell, equal to {—MG\ + NG2) /N g . 
The other reciprocal vector is s = {mG\ — nG2)/N g , 
and we find K = (2n + m)fe/3 + (27V + M)s/3, = 
(n + 2m)fe/3 + (TV + 2M)s/3. 19 In particular, for the 
corrugation along the zigzag direction with A = (6,0), 
both the K and K r points are translated into the origin 
(r-point). 

The numerical calculations for both the flat and the 
corrugated graphene sheets have been performed using 
density functional theory implemented in the ABINIT 



Since any uniform strain respects the inversion sym- 
metry, it can remove the Dirac point degeneracy only 
if it can force the Dirac points to merge, in compliance 
with the Wigner-von Neumann theorem. Merging may 
only take place at a point located at half the reciprocal 
lattice vector, g/2 = (pGi + gG 2 )/2, where p,q are the 
integers 7 . Obviously, the Dirac point can be moved with 
regards to the reciprocal vectors only when the strain 
comprises the shear r. Therefore, below we will distin- 
guish three important cases of the homogeneous strain: 
(i) pure shear strain when eo = 0, en = r, €22 = — r in 
PJ, (ii) uniaxial tension accompanied by Poisson's trans- 
verse contraction, eo = r(l + <j)(1 — <j), €22 = — en<r, with 
a the Poisson coefficient, and (iii) pure uniaxial tension 
(e r, en 2r, e 2 2 0). 

We have performed detailed analysis with account for 
full atomic relaxation of the above cases (i)-(iii) exempli- 
fied by a graphene sheet under a tensile strain e yy with 
varying conditions for a perpendicular strain. In this 
geometry, the Dirac points merge in the case (i), when 
e yy — ~ e xxi an d the strain e yy reaches 23.5%, while in the 
cases (ii) and (iii) for e yy exceeding 26.5 and 27.5%, re- 
spectively. Since the rate of a*— band lowering increases 
in passing from (iii) to (ii) and then from to (ii) to (i), the 
gap is unable to open up in the cases (i),(ii) at all. Al- 
though the gap does open up in the case (iii), its maximal 
value is modest and does not exceed 100 meV since the 
antibonding <r*— band quickly moves down in energy and 
closes the gap (see Fig. [Tl cf. 9 ). The suppressing role of 
<j*— band can be traced back to the fact that the uniaxial 
strain changes the bond lengths and bond angles between 
C-C bonds along the stretching direction and perpendic- 
ular to it. The sp 2 hybrid tends to split into two sp\\ 
hybrids and the perpendicular to the stretching direction 
p± chain states get progressively weaker coupled between 
themselves. As a result, the a*— band derived from p± 
states quickly moves down at the point fc/2^thus fully 
suppressing the gap. 
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FIG. 1: (color online), ((a) The evolution of the gaps be- 
tween the valence ty and the conduction ty* and a* bands 
under uniaxial strain schematically shown in panel (b). The 
Dirac points move along the straight trajectory in BZ before 
merging (c). The bands along the merging direction from Gi 
to G2 are shown in the inset in (a) with the left (right) inset 
showing the bands before (after) the merging of Dirac points. 




FIG. 2: (color online), (a) The evolution of the gaps between 
the valence ty and the conduction ty* and a* bands under 
pure shear strain shown schematically in panel (b). The Dirac 
points move along complex trajectory in BZ before merging 
(c), where f 1 = (Gi - G 2 )/3, f 2 = (2Gi + G 2 )/3. The gap 
between ty* and ty states opens up, but not before the a* — ty 
gap closes (inset in (a)). 



The above results are obtained in the local density ap- 
proximation (LDA) for the exchange-correlation energy. 
Using the generalized gradient approximation (GGA) 
instead does not change the results significantly, even 
though GGA usually favors gap opening. For instance, in 
the case (ii) GGA gives the critical value for the gap open- 
ing € yy = 25.3% with the maximal gap value of about 100 
meV. Keeping this in mind, below we present only the 
LDA results. 

Scanning all possible deformations (i)-(iii) correspond- 
ing to different angles <\> shows that the gap opening is 



most problematic in the case of deformations involving 
stretching along the armchair directions. Here, the D, 
D' points should move the longest distance in fc-space 
to meet each other, and the latter can happen only well 
above the mechanical breaking strain. At the same time, 
the pure shear deformation conjugated to C44 elastic 
modulus {<j) = 7r/4, en = —622) was found to be the 
most favorable for the Dirac points merging, in agree- 
ment with Ref.[T2l Evolution of the band structure under 
such a strain is shown in Fig. [2] In this case, the conduc- 
tion a*— band minimum moves down (inset in Fig. [2^i) 
and intercepts merging of two Dirac points when shear 
strain reaches 17%, while the merging occurs only when 
graphene sheet is strained by about 22%. One concludes 
that the homogeneous strain is unable to open up a gap 
in graphene, contrary to the claim in Ref. [T2l 



IV. GAP IN GRAPHENE SUBJECT TO TWO 
TYPES OF ID CORRUGATION 

Fortunately, the gap can be opened up in all but the 
armchair directions by the corrugation, with the maxi- 
mal gap of about 0.5 eV for the zigzag direction. One can 
distinguish two qualitatively different regimes depending 
on a phase of the corrugation tp: (a) tp = 0, the inver- 
sion (and time reversal) symmetry is preserved, and (b) 
(f ^ 0, the inversion symmetry is generally broken. For 
the case (a), we have calculated corrugation along the 
zigzag direction A = (6, 0), (p = in Fig. [3] Such a per- 
turbation is momentum dependent preserving time rever- 
sal and inversion symmetry and therefore can not auto- 
matically open up the gap in the graphene sheet 25 unless 
the Dirac points are merged at some point g/2. 7 In the 
case under consideration, the Dirac point moves on the 
symmetry V — X line in the folded BZ, where X = fe/2, 
and 'bounces' off the X-point when A/ A is close to 0.12. 
Finally, the Dirac points merge at the T— point. The gap 
opens up when A/ A exceeds a critical value of 0.13 and 
quickly reaches a substantial size of 0.5 eV, as shown in 
Fig. [3] for A/X = 0.15 (see also Fig. [6^). The wider max- 
imal gap relative to the case of a stretched flat lattice is 
simply due to the fact that now the a* band is affected 
significantly less. 

For possible applications, it is important to establish 
if the gap can be opened up when corrugation runs in a 
general chiral direction with an arbitrary phase tp. Band 
structure of graphene subject to a corrugation along a 
'chiral' direction, A = (-2,3), <p = 3tt/19, A/X = 0.14 is 
shown in Fig. [4] Since ip ^ and the inversion symmetry 
is broken, the gap opens up at any i^0. However, the 
gap is very small if the 'pseudo' Dirac points remain well 
inside the BZ, which is the case for small corrugations. 
When the Dirac points approach the BZ boundary with 
increasing corrugation, the minimum of the conduction 
and the maximum of the valence band move apart in 
/c-space. Consequently, the initially small gap grows to 
(indirect) gap of 0.5 eV, Fig. [4] 
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FIG. 3: (color online). Bands in graphene (T — X direction) 
corrugated along the zigzag direction with A = (6,0), ampli- 
tude to period ratio A/X = — 0.15 and the phase (p = 
preserving the inversion symmetry. The Dirac ponts merge 
and the gap opens up when the corrugation is slightly smaller 
than A/X = 0.13. 
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FIG. 4: (color online), (a) Gap in graphene corrugated along 
the chiral direction A = (—2,3), with A/X = 0.14 and the 
phase ip = 37r/19 breaking the inversion symmetry. The Dirac 
cones detach at A ^ and after that their remainders evolve 
differently with corrugation forming the indirect gap (b),(c). 

We should stress that in several aspects our results are 
in contrast with the gap behavior claimed on the basis of 
the continuous flexo-' magnetic' models—* 26 * 27 *. First, the 
latter are oblivious of the presence of the inversion sym- 
metry or its lack thereof 2 ^. As a result, they are missing 
the important increase in the gap size from trivially ex- 
pected minute value in the symmetry-breaking case to 
large values when the remainders of the Dirac cones ap- 
proach the BZ boundaries, Fig. [4] Second, even if the 
inversion symmetry is preserved, the gap can still open 
up when corrugation exceeds some limiting value. Such 
a gap opening does not require any periodic gating cor- 
related with the corrugation suggested iri^l And third, 



it has been claimed that one-dimensional corrugations 
lead to creation of flat bands and a build up of the den- 
sity of states (DOS) at the Dirac (neutral) pointP^ 7 ! 
These partially flat bands are believed to be analogs of 
the zero-energy Landau levels because the ripples affect 
the electrons like an effective magnetic field. Below, we 
show that in reality no flat bands occur since the Dirac 
points are stable up to the deformation where the system 
transforms into a metal. 

As in Ref|26j we restrict ourselves to the limiting case 
of the armchair direction A = (n, n) with the inversion 
symmetry being preserved. Such a corrugation has a 
period of 2n&o, where bo = y/3a/2 (the unit used in 
Ref. [26]). When A/X = 0, there are 2n conduction 7r* 
bands and 2n valence tt bands in the V — X direction per- 
pendicular to A and containing a Dirac point (in the su- 
percell BZ) 17 . Of these 2n bands (conduction or valence), 
two are nondegenerate and n—1 are doubly degenerate^. 
The 7r* and tt bands that form a Dirac cone are always 
nondegenerate (except for the Dirac point itself). When 
the corrugation is induced (A/X ^ 0), all the doubly de- 
generate bands split due to breaking of the mirror sym- 
metry. As A/X is further increased, the lowest in energy 
split 7r* band and the highest in energy split tt band start 
moving toward the Fermi level, due to a change in the 
tt z — tt z hopping matrix elements and increasing <j* — 7r* 
and a — tt hybridization. 

To be more specific, consider the ripples with n = 20 
and A = 40 6o 7 Fig. [5] Here, the closest to the Dirac cone 
are the conduction and valence bands residing inside the 
Dirac cone almost symmetrically (A/X = 0, Fig.[5|. Such 
a relative position of the bands can be easily predicted 
by using the "band-folding" procedure. As the A/X in- 
creases, the degenerate 7r* and tt states split and repel 
each other, so that the lower/upper singly degenerate 
7r* jix state lowers/rises in energy. By approaching the 
Dirac cone, these two states considerably deform the lat- 
ter. At some critical moment (A/X ~0.145), the moving 
7r* and 7r bands reach the Fermi level leading to a peak 
in the DOS at the Fermi level (due to the evident van 
Hove singularities). In Ref. |26l this situation was inter- 
preted as a formation of a flat pseudo-Landau level. In 
reality, however, we face here a peculiar electronic topo- 
logical transition of the type semimetal-metal when the 
bottom of the 7r* band and the top of the tt band pass 
the Fermi level simultaneously. It is clear that such a 
transition is accompanied by electrons flowing from the 
7r to the 7r* band. 

So far, we considered sine- wave deformations with only 
flexural (out-of-plane) atomic displacements. Such de- 
formations are accompanied by the stretching A 2 k 2 /4 = 
(nA/Xy which can be relieved by allowing the atoms to 
relax at a given amplitude A. 'Annealing' increases the 
initial quenched ratio A/X because A shrinks trying to re- 
store the initial nearest-neighbor spacing. Since, on the 
other hand, the pace of movement of the Dirac points 
with A in the relaxed structures is slower, the critical 
corrugation (A/X) c for the gap opening becomes consid- 
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FIG. 5: (color online). Bands in graphene corrugated along 
the armchair direction with A = (20,20), with amplitude to 
period ratio A/X = — 0.15 and the phase ip = preserving 
the inversion symmetry. The k y direction is perpendicular 
to the rippling direction. Note that the bands are defined in 
order of energy. 

erably higher. For example, the annealing of the corru- 
gated structure A = (6, 0) increases the critical corruga- 
tion from 0.12 to ~ 0.5, Fig. [6^i. It is important that 
in the relaxed structures, the a* band is affected signif- 
icantly less than in the case of a stretched lattice and, 
therefore, one may open up a wider gap (~ 1 eV). 

V. DISCUSSION AND CONCLUSIONS 

It is clear from the above that the periodicity of the 
corrugations must not necessarily be commensurate with 
the undistorted graphene lattice or even sine-like for the 
gap opening. The commensurability between the defor- 
mation and the lattice itself would be important if the 
energy gap were due to mixing of electronic states be- 
longing to two different valleys K and K f . However, in 
the case of ID periodic corrugations this does not occur 
even if the corrugations provide a momentum transfer of 
k = K — K' (as in Fig.l) The gap opens only due to 
(i) breaking of the inversion symmetry and/or (ii) merg- 
ing of two inequivalent Dirac points, as explained above. 
We simulated, for example, the situation when a corru- 
gated graphene is placed on a hexagonal boron nitride 
(h-BN) and only a fraction of its atomic lines are in reg- 
istry with h-BN. Such a partially coherent connection 
with the substrate maintains Gaussian-like periodic cor- 
rugations that we may call a 'mesh' (see the insert in 
Fig. [6)3). In Fig. J6Jd, we present the results for a partic- 
ular case when the graphene is deformed in such a way 
that its initial period A = (6, 0) becomes commensurate 
with the BN period A = (4, 0). We see that the corruga- 
tion leads to a gap of 0.20 eV. In practice, the periodic 
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FIG. 6: (color online). Band gaps in corrugated periodic 
graphene sheets, (a) For quenched and annealed sinusoidal 
corrugations with A = (6, 0) as a function of the corrugation 
parameter A/X. The critical corrugation (A/X) c significantly 
increases in passing from the quenched to annealed ripple 
structures. The annealed structure at (A/X) c is shown in the 
insert, (b) For a corrugated graphene placed on a hexagonal 
boron nitride (h-BN) substrate. It is supposed that the con- 
nection with the substrate is only partially coherent: within 
the period only four graphene atomic lines are in registry with 
h-BN lattice (shown by vertical dotted lines in the insert). 



deformations can be achieved in various ways. One can 
imagine a situation when the edges are not clamped and 
a graphene membrane bends in accordion-like fashion. 
This situation corresponds to the 'annealed' corrugations 
where the initial nearest-neighbor lengths are mainly pre- 
served. The 'quenched' corrugations can be induced by 
depositing graphene on a substrate like Si02 28 or the vic- 
inal surfaces with regularly spaced steps, like Au(788)P^. 
In the case of Si02, for example, the interaction energy 
between the graphene sheet and the substrate is suffi- 
cient to overcome the elastic energy needed for graphene 
to conform to the Si02 surface profile 28 . 

In conclusion, we have shown that any practical homo- 
geneous deformation cannot open up a gap in graphene 
sheets. At the same time, the inhomogeneous deforma- 
tion can open up a significant gap (~ 1 eV) rather in- 
dependently of direction and form of a corrugation with 
an exception of an armchair direction. The present gap 
opening does not require any periodic gating correlated 
with the corrugation^ either. It is worth noting that in 
the limit of a clean undoped graphene sheet, the Coulomb 
interaction between electron and hole carriers may be- 
come significant if the Coulomb coupling constant ex- 
ceeds some critical value on the order of unity and may 
facilitate an excitonic gap in the spectrum^]. However, in 
practical cases the Coulomb interaction may be screened 
by (unintentional) doping. Since a multilayer graphene 
( ML G) /graphite is a stack of weakly coupled graphene 



6 



sheetd^, most of above results, therefore, should apply to the inhomogeneously deformed MLG as well. 
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